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Abstract 

We consider a one-dimensional continuum gas of pointlike positive and negative 
unit charges interacting via a logarithmic potential. The mapping onto a two- 
dimensional boundary sine-Gordon field theory with zero bulk mass provides the 
full thermodynamics (density-fugacity relationship, specific heat, etc.) of the log- 
gas in the whole stability range of inverse temperatures (3 < 1. An exact formula 
for the excess chemical potential of a "foreign" particle of an arbitrary charge, put 
into the log-gas, is derived. The results are checked by a small-/? expansion and at 
the collapse /3 = 1 point. The possibility to go beyond the collapse temperature is 
discussed. 
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1 Introduction 

The model under consideration is a two-dimensional (2D) two-component (TC) classical 
Coulomb gas confined to a ID manifold, the circle of radius i? or its i? — » oo limit - the 
straight line. It is usually named as the symmetric ID TC log-gas. The model consists 
of a mixture of mobile pointlike particles {j} of positive and negative, say unit, charges 
{qj = ±}, localized at continuous angle positions {^Pj} on the circle of length L = 2nR. 
The interaction energy reads 

E{{qj, <^j}) = IjIkviVj, Vk) (1.1) 
j<k 

where v{ip,ip') is the 2D Coulomb potential with imposed periodic boundary conditions, 

t;(^,V50 = -ln{(^)|e'^-e^^'|} (1.2) 

The length constant tq, which fixes the zero point of energy, will be set for simplicity to 
unity. In the thermodynamic limit R —>■ oo, which is of interest in this paper, introducing 
the straight line position variable x G (—00,00) via x = (fR, the interaction potential 
( p,.2D takes the familiar logarithmic form 

v{x,x') = — \n{\x — x'\) (1.3) 

For pointlike particles, the interaction Boltzmann factor of a positive-negative pair 
of charges at distance x, x~^ where (3 is the inverse temperature, is integrable on a 
ID manifold at small x if and only if /5 < 1 (at large x, there is no problem because 
the interaction is screened by the system), with (3 = 1 being the collapse point. In 
its conductive phase, the model exhibits poor screening properties, a typical feature of 
a D-dimensional Coulomb system confined to a domain of dimension D — 1 ll|]. As a 
consequence, the leading non-oscillatory term of the large-distance asymptotic decay of 
the charge-charge correlation is algebraically slow ~ — l/[/3(7ra;)^] 0. 
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The ID TC log-gas, with ± charges required to ahernate in space, is equivalent to the 
impurity Kondo problem [§] - 0. 

The ID TC log-gas, without any restriction on the charge order, is related to the 
problem of quantum Brownian motion of a particle in a ID periodic potential [^] - ||^. 
The model was also studied in connection with the problem of non-equilibrium quantum 
transport through a point contact in a ID Luttinger liquid ^ - [11|, having a realization 



in resonant tunneling-transport experiments between edge states in fractional quantum 



Hall effect devices fi^ . 



The ID TC log-gas with a hard core, or with some other short-distance regulariza- 
tion of the Coulomb interaction potential, undergoes a ID counterpart of the Kosterlitz- 
Thouless transition from a conducting phase to a low-temperature insulating phase around 
j3 = 2. An approximate analytic study of the gas with log(l + interaction was given 



in ref. [0|. The lattice version of the model was exactly solved at /?= 1,2,4 and its 



conductor-insulator phase diagram was conjectured in refs. Ijl^ - |T6| . 

Recently |T^, the bulk thermodynamic properties (density-fugacity relationship, free 



energy, internal energy, specific heat, etc.) of the infinite 2D TC plasma have been 
obtained exactly, in the whole model's stability range of temperatures, by mapping the 
plasma onto a 2D bulk sine-Gordon field theory and then using recent results about that 
integrable field theory. The aim of this paper is to derive the full thermodynamics of the 
confined plasma, the ID TC log-gas, via its relationship to the 2D boundary sine-Gordon 
field theory with zero bulk mass. The particle nature of the statistical model permits one 
to check the obtained behaviour of thermodynamic quantities close to the collapse f3 = 1 
point, and to suggest a possible analytic extension of the results beyond that collapse 
point. 

The mapping onto the 2D boundary sine-Gordon model is outlined within the grand- 



canonical formalism in section 2. To give to the mapping a precise meaning, we analyse 
the short-distance behaviour of the pair distribution function for two basic combinations 
of particle charges. The relationship between the ordinary statistical mechanics |TE], 
T9|] and the Conformal Perturbation theory [^] is established. As a byproduct, for an 



"external" (or, probably more adequately, "foreign") particle with an arbitrary charge 
put into the log-gas, a relation between its excess chemical potential and the one-point 
expectations of the exponential boundary field in the sine-Gordon model is derived. 
Section 3, devoted to the derivation of full thermodynamics of the ID TC log-gas, is 



based on an exact fromula for the above one-point expectations ^T[, obtained by using a 
"reflection" relationship between the 2D boundary Liouville and sinh-Gordon theories. 

Section 4 brings a discussion about a possible analytic extension of the acquired results 
to the collapse region 1 < (3 < 2. In this region, one could attach to the particles a small 
hard core a in order to prevent the collapse, then calculate thermodynamic quantities, and 
at the end take the a — >■ limit. In this limit, while the free energy and the internal energy 
per particle diverge, the specific heat, truncated correlation functions, etc., are expected 



to remain finite An analytic continuation of the formula for the specific heat (at 

constant volume) predicts an infinite sequence of phase transitions from the conducting 
plasma phase /5 < 1 to the insulator region (3 >2. In two dimensions, such a phenomenon 



was predicted in ref. but later denied in ref. |24 



Whenever possible, the results are checked by a systematic small-/? (high-temperature) 



expansion, using a renormalized Mayer expansion technique developed in refs. ||T^, |25 



26|1 (see Appendix), and close to the collapse [3 = 1 point, by applying arguments in the 



spirit of an independent-pair approximation p2| , which assumes a dominant contribution 
from almost collapsed positive-negative pairs of charges to the configurational integral. 
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2 Mapping onto the boundary sine-Gordon theory 



Let us first examine a general TC plasma of g = ±1 charges in an infinite 2D space of 
points r = (x, y). We will work in the grand canonical ensemble, with position-dependent 
fugacities z+{r) and Z-{r) of the positive and negative particles, respectively. In infinite 
space, — A/(27r) is the inverse operator of the Coulomb potential — ln(|r|). Using the 
standard procedure (see, e.g., ref. |2^), the grand partition function S of the plasma at 
inverse temperature (3, considered as the functional of particle fugacities, can be turned 
into 

JV(f)exp [/dV ff 0A0 + z+(r)e'v^'^ + z_(r)e-'v^'^) 

^ (2.1) 



^ Z:, z_ 



JV<Pexp (/dV^. 

Here, 0(r) is a real scalar field, / T>(f) denotes the functional integration over this field and 
the fugacities are renormalized by a self-energy term. The consideration of 



( ^) = Zgix)5{y) 



(2.2) 



confines the charges to the x-line. 



[z+,z. 



JVcpexp < 


/dV(^0A0) 


+ /f^dx 


2+(a;)eV 




} 


/P0exp (/d2r^0A0) 





(2.3) 



where 0b(x) = (f){x,y = 0) is the boundary field. 

To reformulate the field theory (p.3[) as a boundary problem, in formal analogy with 



ref. 281, one introduces two new fields 



h{x,y) 
^o{x,y) 



-j= [0(x,t/) + 0(x, -y)] 



V2 



{x,y) -0(x, -y)] 



(2.4a) 
(2.4b) 



defined only in the upper half-plane y > 0. The even field has a Neumann boundary 
condition dy(j)e\y=o = and the odd field has a Dirichlet boundary condition (f)o\y=o = 0. 
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It holds 

/ dV0A0 = / dV(0eA0e + 0oA0o) (2.5) 

J Jy>0 

The odd field, contributing only by its free-field part 0oA0o, disapears from ( |2.3| ) by 

numerator-denominator cancelation. By integration per partes, the term ^e^^e can be 

rewritten as — (V0e)^; with a vanishing contribution from the boundary. Considering that 

4'b{x) = 4>e{x^y = 0)/\/2, and renaming then 0e as 0, ( p.3|) transforms to 

^ . _ /P0exp(-g,G[^+,2_]) 

/I)0exp(-5.G[O,O]) ^'-'^ 

with the action 

/OO /"OO 1 /"OO p -| 

dx dy—{V<pf~ dx z+(a;)e^^'^^ + z_(x)e-^''^^ (2.7a) 
-OO JO 47r J-oo L J 

/? = 262 (2.7b) 

where again (I)b{x) = (t>{xiy = 0) and the (/)-field has the Neumann boundary condition. 
For uniform and equivalent charge fugacities, z+{x) = z^{x) = z, one gets the boundary 
sine-Gordon model with zero bulk mass, 

dx dy— (V0)2-2z/ da:cos(60B) (2.8) 

-OO JO 47r J-oo 

The sine-Gordon representation of the multi-particle densities can be obtained from 
the functional generator ( p.6|) , (|2.7|) in a straightforward way. The density of particles of 



one sign is 

51n; 



= z,{e'^'"f'^)sG (2.9) 

Zq(x) = Z 



5Zq{x) 

where {■■■)sg denotes the averaging over the sine-Gordon action (^.8]). Here, although 



the charge symmetry is considered, i.e., n+ = n_. = n/2 (n is the total number density of 
particles), we leave the g-subscript in order to make transparent identities like {e*'^'^)sG = 
(e~''"^^)sG- This identity is a special case of a general symmetry relation 

{e''"^B)sG = (e-'°'^^),G, a arbitrary (2.10) 



which results from the invariance of the sine-Gordon action ( p.8| ) with respect to the 
transformation —(f). The pair distribution function gg^q>{x,x') is given by 

6^E 



9q,q'{\x-x'\) 







nj 


[Ug, 




(Zg^ 


rig) 


\ng, 



Ugi J E 6zg{x)6zq'{x') 



Zq{x)=Z 



^\qb<t>B(x) iq'b4>B{x') 



)sG (2.11) 



and so on. 

In statistical mechanics, the short- distance behaviour of the pair distribution function 
9q,q'{\^ — x'l) is dominated by the Boltzmann factor of the pair potential, 

~ Cq^q'\x - x'\'^'"'' as \x - x'l ^ (2.12) 

(provided (3 is small enough; see below). The prefactors Cq^qi are related to a free energy 
difference [|T^, [jl9[. For q' = —q, one has 



Cq,.q = exp [(3 (/i^^f + /i!.^)] (2.13) 

where fi^^ is the excess chemical potential of species q. Strictly speaking, relations (|2.12|) 
and ( |2.13| ) are valid only for /9 < 1; at /? = 1, due to the collapse of positive- negative 



pairs of charges, /if" oo. On the other hand, for the truncated two-body densities 
ng j,{x, x') = nqnqi[gg^gi{x, x') — 1], the short-distance asymptotics analogous to that given 
by eq. ( p.l2| ) with q' = —q takes place also for /3 > 1. For q' = q, one has 

C,,, = exp [/? (2/xf - (2.14) 

Here, we use an extended definition of the excess chemical potential for "foreign" ions 
with arbitrary charges, put into the underlying electrolyte: (i'q = reversible work which 
has to be done in order to bring a foreign particle of charge Q from infinity into the 
bulk interior of the ID TC plasma of unit ± charges (the consequent breaking of the 
system neutrality has a negligible effect in the thermodynamic limit). This quantity is 
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of evident importance in chemistry. The interaction Boltzmann factor of the Q charge 
and an opposite unit charge at distance x, is integrable at small x if and only if 

\Q\[3 < 1. The stability region for /Xq therefore is expected to be /5 < l/IQI with /Xg 
diverging just at /? = l/IQI- As a consequence, relations ( p.l2| ) and ( |2.14| ) are valid only 



for (3 < 1/2. For (3 > 1/2, the two-point correlator oc |x — x'\^^^; the analysis of this 
change in the short- distance expansion goes beyond the scope of the present work. With 
regard to the thermodynamic relation 



In^i =/3/i- (2.15) 



we conclude that, for g = ±, it holds 

{e'^'"f-),G = exp{-f3f,'^^) (2.16) 

and 

(ei9fe0BWe-igb<^s(x'')^^^ _ |a;_a;'|-26^ as|x-x'|^0 (2.17a) 

^^iqb<PBix)^iqb<PBix')^^^ _ g'^^^^ | X - xf ''^ as|x-x'|^0 (2.17b) 

where we have used the equality P = 26^, see formula ( [^.7b| ). 

In quantum field theory, the sine-Gordon model ( p.8|) can be regarded as a conformal 
field theory perturbed by the boundary cos-field. The short distance expansions for multi- 
point correlation functions are then obtainable by using the Operator Product expansions 
(OPE) The vacuum one-point expectations, in our case (exp(ia0B))sG with an 

arbitrary value of a, are the basic objects of the OPE scheme which contain the whole 
nonperturbative information about the system. For the product of two primary fields 
e'"i'^(xi)e'''2</.(^a,2), the OPE has the form 

oo 

e''^i<^«(xi)e^'^^^^(x2) = [Q°a2(ki-a;2|)e'("^+"^+"')^^(x2) + ---] (2.18) 



where the coefficients C are computable within the Conformal Perturbation theory |20[ . 



By successive apphcation of (|2.18| ) the short-distance behaviour of any multi-point cor- 



relation function of the exponential field can be reduced to one-point functions. For the 
present model, the leading short- distance term of the two-point function reads 

^giqfe</.sWeVfe0s(x'))^^ ~ {S"^'''^^^^) sg\x - x'\'"''''^' as \x - x'| ^ (2.19) 



For q' = —q, formula ( p.l7a| ) is reproduced. For q' = q, the combination of (|2.17b|) and 
( p.l9| ) implies the relationship 



exp 



{-Pf^i) = {e'"'"^^)sG (2.20) 



The above formalism generalizes straightforwardly to Q-particle distribution functions 
{Q positive integer), with the result 

exp = (e^«^^-).G (2.21) 



Eqs. ( p.l6| ) and ( p.20|) are the lowest Q = 1 and Q = 2 cases, respectively, of this formula. 



We expect the validity of relation (p.21|) also for noninteger values of Q. An analogous 



formula can be derived for the infinite (unconfined) 2D TC plasma. 

3 Thermodynamics 

The 2D boundary sine-Gordon theory with the action ([2.8|) has a discrete symmetry 



0^0 + 27Tn/b with any integer n. In the domain < 6^ < 1 this symmetry is sponta- 
neously broken and the theory has infinitely many ground states |0„), characterized by 
the associate expectation values of the field 0, (0)„ = 2Txn/h. One has to choose one 
of these ground states, say |Oo), and consider (■ ■ ■)sg = (■ ■ ■)o- The parameter z, which 
renormalizes multiplicatively, gets a precise meaning when one fixes the normalization of 
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the field cos(60b). The conformal normahzation ||2T[, based on the OPE scheme (|2.18|) 
and consistent with formula (|2.17a| ) derived via the ID TC log-gas, corresponds to the 
short- distance limit of the two-point function 



)sG ^ \X -X 



/|-2a2 



clS I X X I ^ 



(3.1) 



Under normalization (|3.1|) , the expectation value of the exponential boundary field was 



obtained in ref. [pT 





aV(l-62) 1 


/■°° dt 


_r(62)_ 


exp < 


Jo T 



(e* - 1 + e*(i-^') + e-^'*) sinh'' (abt) 
2 sinh(62t) sinh(t) sinh((l - b'^)t) 

o / 1 



—a 



+ e" 



(3.2) 



^sinh((l - 62)t) 

where |Re 2a| < 1/6. This result was derived under assumption of a "reflection" relation- 
ship between the 2D boundary Liouville and sinh-Gordon theories (with zero bulk mass), 
and the consequent analytic continuation of the latter theory in the 6-parameter to the 
boundary sine-Gordon model. For the case of special interest a = b, using the integral 
representation of the Gamma function 



In r(x) 



°o dt 

— e 
t 



X-1 + 



-{x-l)t _ I 

1 -e-* 



Re x > 



(3.3) 



after some algebra, relation ( p.2|) takes a simpler form 

1 „/l-262\_/ 62 



sG 



.r|l:i£,r 



47r3/2(l -62)^ y2-2by \2-2b\ 
On account of the symmetry relation ( |2.10| ), it holds 



2nz 

W) 



1/(1-62) 



(3.4) 



(e'''<^-).G = ^|- lim ^InS^ 

2 OZ L-*oo L 



where L is the length of the system. Since Siiz = 0) = 1, one has 



(3.5) 



ir^ z = 2^^ ^ (2^ 



27rz 
f(62j 



l/(l-fe2) 



(3.6) 
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Relating the sine-Gordon parameters with those of the TC log-gas, (3 = 26^ from 
( p.7b|) and n/{2z) = (exp(i60B))sG according to ( ^^ , formulae in the above paragraph 
imply the explicit n — z relation for the plasma: 



n 



1-/3/2 



27r 



r(/5/2) 



1 1-/3/2 



(3.7) 



The small-/? expansion of the rhs of 

^1-/3/2 



reads 



z 



2/3^/^exp^(C + lnvr)| + l[.r^ + 7C(3)] (| 



+ 12 



TT^ + 6C(3) 



1^2 



0(/5^ 



(3. 



where C is the Euler's constant and C is the Riemann's zeta function. Formula ( |3.8D is 
checked in the Appendix by using a renormalized Mayer expansion in density. For fixed z^ 
relation ( p^) tells us that the particle density n exhibits the expected collapse singularity 
as /? — >• 1~: 

(3.9) 



n 



1-/3 



The form of this singularity can be reproduced indirectly by using a perfect-screening sum 
rule for the one-body densities, 



da; 



(3.10) 



where the truncated two-body densities rv L{x,x'^ = nqnqi[gq^qi{x,x') — 1]. For (3 ^ V 



the integral in (|3.1CI|) is dominated by the short-distance behaviour of nq^_q{x) ~ z'^\x\ ^ 
[see relations (|21^ ), (|27[3| ) and (pJ[5|)]. Then, 



n 
2 



dx 

-A Ixl^ 



2z^ 
l-(3 



1-/3 



2^2 



1-/3 



+ 0(1) 



(3.11) 



where A is a length over which the Coulomb interaction is screened by the system, and 
the singular behaviour (|3.9|) is justified. 
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To get the full thermodynamics of the ID TC log-gas, we pass from the grandcanonical 
to the canonical ensemble via the Legendre transformation 



FL{|3,N) = nL + ^lN 



Here, the total particle number = riL, the chemical potential 



(3.12) 



(3.13) 



and the grand potential Vt is defined by 



j3nL = \nEL = L\l--\n 



(3.14) 

where we have combined eqs. ( p. 41) and ( p.6|) . Note that the thermodynamic relation 
InS^, = 13 PL (P is the pressure) implies the well-known equation of state 



(3P 



(3.15) 



The dimensionless specific free energy f,f = PFl/N, is then written as 



(-C 



-1 + Inn + In M - ^ 



2 2 I 2 



inrfin^Umr^^/^ 



(3.16) 



According to the elementary thermodynamics, the (excess) internal energy per particle, 
u^^ = (E), and the (excess) specific heat at constant volume per particle, = Cy/N, 
are given by 



u 



(3.17a) 
(3.17b) 



For Cy , one gets explicitly 



■-v 



4(2-/3)3 
-2(/5-2)2 + 2^ 



2^' { ^1 2)¥ ( ^ 



/5/2 



2-/5. 



(3.18) 
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where ip{x) = d[lnr(x)]/dx is the psi function. The series representation of its first 
derivative reads 



1 



(3.19) 



The expansion of Cy/kB around the collapse point /? — 1~ results in the Laurent series 



„ex 



1 



kB 2(1-/5)2 2(1-/5) ^ 1,2 ^ 24 



+ 0(1-/5) 



(3.20) 



The leading term of this series can be reproduced by regarding that, under assumption 

/(/5, n) ~ i ln(l - /5) - 1(1 - /5) In A as /5 1" (3.21) 

and then applying ( p.lTbj ). 

As concerns the excess chemical potential of a foreign particle of charge Q in the 



plasma, setting a = bQ and 6 = /5/2 in ( p. 21) , relation (|2.21|) gives 



ex 

Q 



2-/5 



In 



2/3/2 



r(/5/2) 



+ 



°° dt 

T 



(e* - 1 + 



,t(l-/3/2) 



+ e-^*/2) sinh2(g/5t/2) 



2sinh(/5t/2) sinh(t) sinh((l - /5/2)t) 
/5Q2 / 1 



+ e- 



(3.22) 



2 Vsinh((l -/5/2)t) 

In the limit of large t, the first integrand behaves like (1/t) exp[(/5|Q| — l)t], and therefore 
the integral is finite if /5 < 1/|Q|. As was argued in the previous section, this is the true 
stabihty region for /Xq. For \Q\ < 1, /ig remains finite also in a subspace of the collapse 
region, namely in the interval 1 < (3 < 1/\Q\- Everything what was said above holds only 
when (3 < 2. Approaching the conductor-insulator phase transition at point /5 = 2, the 
first term on the rhs of (|3.22|) diverges for an arbitrary Q. 



4 Discussion 



All results in this paper were derived under the assumption that the system is stable 
against collapse, i.e., in the range of inverse temperatures j3 < 1. The only exception 
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is represented by the excess chemical potential /Xq of a "foreign" Q-charged particle 
embedded in the plasma: if \Q\ < 1, relation ( p.22|) applies also to the collapse region, 
up to /3 = 1/|Q|. A natural question arises whether there is a possibility to make an 
analytic continuation of the explicit result for some bulk quantity into the collapse region 
1 < P < 2. The best candidate for such a continuation is the specific heat, formula ( p.l8|) . 



which might be in the sense explained in the Introduction a well defined finite quantity 
also in the collapse region. 



Before going further, let us recall the work of Gallavotti and Nicolo ||2^ concerning the 
infinite 2D TC Coulomb gas of unit ± charges. In two dimensions, the particle collapse 
occurs aX j3 = 2 and the Kosterlitz-Thouless phase transition takes place around (3 = 4:. 
By studying the Mayer series of the specific grand potential in fugacity it was proven in 



ref. |2^ that each term of the series converges in the insulator region (3 > A. For /3 < 4, 



the existence of infinitely many thresholds 



2D: /?^ = 4(1--1-), iV = l,2,... (4.1) 



2N. 

was observed: only the Mayer series' coefficients up to order 2N are finite if /? > j3i\f. Points 
{Pn} were conjectured to correspond to a sequence of transitions from the pure multipole 
insulating phase (/3 > 4) to the conducting plasma phase {(3 < 2) via an infinite number 
of intermediate phases. Although the mathematics used in p3l was quite complicated. 



there exists a simple explanation of the above findings |31|. For a neutral choice of 
positive and A^ negative charges in a disk of radius R, the most relevant contribution to 
the coefficient of the z"^^ term comes from the complete-star integral, after approximating 
the Mayer function by exp(— 

1 rR ^ Y\ ■ \r- — r-\^\r^. — 

ndv.dv:iii^4^^^^i^^ (4.2) 



i=l I ' J\ 

where cr is a small hard-core diameter and r (/) denote spatial coordinates of positive 
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(negative) charges. Rescaling Rsi and Rs[, one gets (4.2) = ^ 

f{cr/R) which diverges as i? ^ oo just when f3 < (3^- In two dimensions, the formula for 
dy/kB in the region /5 < 2 is presented in ref. [FZl, relation (56). The only source of the 



singularities of Cy/ks, extended into the region 2 < /? < 4, is the term oc sin ^(7r/?/(4— /?)), 
which gives a double pole at 



/^^ = Ml-]^) N = l,2,... (4.3) 



For = 1, 3, 5, . . ., these singular points coincide with the ones in ( [4.1[ ). However, there 
are additional divergencies of Cy/fc^ for = 2,4,6,... when Cy/Zc^ — oo, which is 
an unacceptable thermodynamic behaviour. Therefore, in 2D, the analytic extension of 
the formula for d^/ks is meaningless, what supports the arguments of Fisher et al 
indicating the total absence of any intermediate phase at nonzero particle density. 

The situation is different in ID. The dominant configuration integral of the ^-series at 
the 2iVth order 

diverges for the line length L ^ oo when 13 < (3^ with 

ID: ^^ = 2{l-^y iV=l,2,... (4.5) 



The singularities of dy/kB-, given by formula (|3.18| ), originate in the region 1 < (3 < 2 



exclusively from the term ip' {{1 — (3) / {2 — [3)) . According to (|3.19|) , ip'^x) has second-order 



poles at X = 1 — N {N = 1,2,...). The corresponding /J^r coincide just with P^, so the 
analytic continuation of (|3.18| ) into the collapse region 1 < P < 2 reproduces exactly the 



singularities suggested by Gallavotti and Nicolo. Moreover, Cy-fkB is always positive in 
the underlying region 1 < /? < 2. It is therefore tempting to conjecture that the multipole 
intermediate phases might exist in ID. 
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The last remark concerns the generahzation of the relation (|3.11|) to the case of a small 
hard core cr: 



n 



- = 2 dx 



r^Q (^"' - ^"') (4-6) 



2 Ja \x\l^ 1-/3 
with A > (7. Then, around [3=1, the free energy 

/^llnf \ (4.7) 

is well defined for /? ^ 1^ as well as /3 — 1+. Applying ( p.l7b| ) and taking the cr — 
limit, one finds 

7^ ~ on ^ m2 both /5 ^ 1- and /5 ^ 1+ (4.8) 

i.e., the leading singular term of the expansion (^.201) admits an analytic continuation 



from /3 < 1 to the (3 > 1 region. This fact supports the above conjecture. 

We hope to motivate numerical simulations of the model under consideration. 

Finding an integrable TC Coulomb gas with some short-range regularization of the 
interaction Coulomb potential might be a realistic task in one dimension. 



Appendix 

We consider the ID TC log-gas confined to a straight line of size L ^ oo. A pair of 
(±)-charged particles z,j interacts via 

v{xi,qi\xj,qj) = qiqjv{xi,Xj) (A.l) 

where the Coulomb potential v{x,x') is defined by (|1.3|). The supposed equality of the 
species fugacities z+ = = z corresponds to a neutral system with homogeneous particle 
densities = n_ = n/2. 
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The renormalized Mayer expansion in density (for details, see refs. [|T^, ||25[, p6| ) 

is based on the expansion of each Mayer function in the inverse temperature /3, and on 

the consequent series ehmination of two-coordinated field circles between every couple 

of three- or more-coordinated field circles; by coordination of a circle we mean its bond 

coordination. The renormalized bonds are given by 
K -I3v 

(A.2) 

where, besides the integration over spatial coordinate x of a field (black) ng(x)-circle, the 
summation over charge g-states at this vertex is assumed as well. For the interaction 
under consideration ([A.l|) , the renormalized bonds exhibit the same charge- dependence, 

K{x,q\x\(^) = qq'K{x,x') (A. 3a) 

where K{x, x') = K{\x ~ x'\) is given by 

Kix) = r ak^p^ (A.3b) 
^ ' 2J-00 \k\ + 7m(3 ^ ' 

By a simple rescaling of integration variable /c, K is shown to exhibit a special scaling 
form 

K{x) = -(3K{7ml3\x\) (A.4a) 



The small- a: expansion of K reads 

K{x) ~ (3 



ln{7m(3\x\) + C - ^n'^n(3\x\ + ... 



(A.5) 



where C is the Euler's constant. At asymptotically large-x distances, by using twice 
integration per partes it can be shown that K decays algebraically as follows 

\j:npxY 
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which exhibits the poor screening properties of log-gases in ID [^,0. 

The procedure of bond-renormahzation transforms the ordinary Mayer representation 
of the dimensionless (minus) excess Helmholtz free energy, denoted as A [n] , into 

oo 

A[n] = ^ ^ + D^'\nl (A.7a) 

s=l 

where 

D(o)= \ + I ! + ... (A.7b) 

is the sum of all unrenormalized ring diagrams and {Z^^'')} represents the set of all remain- 
ing completely renormalized graphs; the first few D*^^)-graphs are drawn in the sketch (11) 
of ref. ||T^. The first term on the rhs of ( |A.7a| ) is fixed to zero by c harge neutrality. The 



second term 1)^°'', eq. ( |A.7b|) , is expressible as follows [|T^ 



= _ / dn' lim \K{x, v!) + ^v{x)\ (A.8) 

2 JO 2:^0 

With respect to ( [A.5| ), one gets 

D'^^^ 3 3n 

— - = ^(nlnn-n) + ^[C + ln(7r/?)] (A.9) 

As concerns the completely renormalized graphs, due to the ± charge symmetry, only 
those Z)*^*) are nonzero whose all vertices have an even bond coordination. The scaling form 
of the dressed bond formula ( |A.4D , permits us to perform the n- and /3-classification 
of a nonzero diagram D'^^\ composed of Ng vertices (each bringing the factor n) and Lg 
bonds. Each dressed bond contributes by the factor —[3 and enforces the substitution 
x' = xnnf] which manifests itself as the factor l/{7m(3) for each field-circle integration 
/ dx. Since there are {Ns — 1) independent field-circle integrations in D^^\ one concludes 
that 

±1_ = n/3^»-^»+M, (A. 10a) 



where is the number 



= 1,(3=1) 



The first nonzero diagram from the sketch (11) of ref. |T7 



IS 




It contributes to the (3^ order, with 



do 



In the next order, only diagram 



1 p dx -4 




survives, and 
(if 



1 



3!(2!)3y_oo vr J-oo vr 
etc. To evaluate d2 and dg) we Fourier-transform K'^{x): 

°° dx 



G{k) 



-oo y 27r 

(l + |fc|) 
|A;|(2 + |A;|) 



ln(l + |A;|) 



In terms of G{k), 



d2 
de 



2!4!7r .y-oo 
1 



dk G\k) 



2tt / dA; G%k) 



3!(2!)%2 

With the aid of the substitution 1 + A; = exp(t) {k and t positive), 

1 1 



d2 

da 



12 23 

1 1 

12 24 



Tx' + 7C(3) 
vr^ + 6C(3) 
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where ( is the Riemann's zeta function. 
Finally, for g = ±, 



ln(^ 

z 



' , 22; ) dn 



A(n) 



Using results of the previous paragraph, one ends up with 



n 



2/3^/2 I (c- + In tt) I + 



exp \{C + Invr) 



+ 



12 



TT^ + 6C(3) 



'/5 



2 12 

, 4 



+ 7C(3) 



(A.16) 



(A. 17) 



in agreement with (|3.8|). 
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